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Abstract: Using the low energy limit of type IIB superstring theory we obtain the non- 
extremal limit of deformed conifold geometry which is dual to the IR limit of large N thermal 
QCD. At low temperatures, the extremal geometry without black hole is favored while at 
high temperatures, the field theory is described by non-extremal black hole geometry. We 
compute the ten dimensional on shell action for extremal and non-extremal geometries 
and demonstrate that at a critical temperature T c there is a first order confinement to 
deconfinement phase transition. We compute T c as a function of 'tHooft coupling and 
study the thermodynamics of the dual gauge theory by evaluating the free energy and 
entropy of the ten dimensional geometry. We find agreement with the conformal limit 
while thermodynamics of non-conformal strongly coupled gauge theories is explored using 
the black hole geometries in non-AdS space. 
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1. Introduction 

At extreme temperatures nuclear matter is best described by a weakly interacting gas of 
gluons and quarks due to asymptotic freedom ]l], gj. The weak coupling allows one to use 
perturbative techniques to study the thermodynamics of the system and at zeroth order in 
perturbation theory, it is best described as a gas of free particles. But the same asymptotic 
freedom implies that the perturbative analysis must break down with decreasing temper- 
ature. As temperature is lowered, the couplings get stronger and color degrees of freedom 
are confined. Thus nuclear matter undergoes a phase transition- from low temperature con- 
fined phase of color neutral constituents to high temperature deconfined phase of quarks 
and gluons. However to analyze the non-perturbative confinement mechanism, one has to 
study the theory on the lattice |4], ||] or resort to effective field theories - both of which 
are successful along with their inherent limitations 1 . 

On the other hand 'tHooft showed that when number of colors N goes to infinity, 



only planar diagrams contribute to the amplitude |10| , 11]. Thus in the large N limit, 
the theory drastically simplifies and one may hope that the thermodynamics could be 
analytically tractable. Furthermore, the structure of the planar diagrams suggest that 



1 For analysis of high temperature phase and phase transitions please consult For a review of 

recent developments in lattice QCD, please consult and for the large N limit. 
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four dimensional gauge theory maybe interpreted as a string theory [12]. This connection 
becomes clear with the realization that gauge theories naturally arise from excitations of 
strings ending on branes [35], while on the other hand gravitons arise from low energy 
excitations of closed strings. By studying the interactions between open and closed strings, 
the correspondence between gauge theory and gravity can be realized. The best studied 
example is the AdS/CFT correspondence proposed by Maldacena |3(| which essentially 
maps maximally super symmetric J\f = 4 SU(N) conformal gauge theory to AdS§ x S* 5 
geometry. In the large N limit, the gauge theory has large 'tHooft coupling N — > oo, g s — > 
0,g s N » 1 and AdS$ x S 5 geometry with small curvature ~ l/L 2 ,L = (g s A r ) 1//4 v / o7 » 1 
has a classical description. Again the appearance of large N has allowed an exact analytic 
description of strongly coupled quantum gauge theory in terms of weakly coupled dual 
classical gravity. 

In [37] an exact proposal was made to compute correlation functions of four dimen- 
sional J\f = 4 conformal SU(N) gauge theory using partition function of supergravity on 
AdS space. To study the thermodynamic properties of the gauge theory, one studies the 
thermodynamics of the dual geometry and for the case of CFTs, this amounts to the study 
of black holes in AdS space. The AdS/CFT correspondence has been quite successful in 
this regard- the scaling of thermodynamic state functions (free energy, pressure, entropy) 
of CFTs with respect to temperature matches with that of black holes [^, [HJ. However 
there are no phase transitions in conformal field theories and the strong coupling regime 
of QCD matter is far from conformal. Thus the study of confinement to deconfinement 
phase transition in QCD requires careful extension of the gauge/gravity correspondence to 
incorporate the conformal anomalies. 

In principal this can be done by placing D branes in ten dimensional geometries and 
then studying the warped geometry sourced by the fluxes and scalar fields arising from 
the branes. Since QCD is non supersymmetric and non-conformal, the first objective is 
to find gauge theories with RG flows arising from D brane configurations with minimal 
SUSY. There have been a lot of progress in that direction: In pl| , [42"[ RG flows that 
connected conformal fixed points at IR and UV was incorporated, and [43] connected the 
UV N = 4 conformal fixed point to a J\f = 1 confining theory. But the model with QCD 
like logarithmic running of the coupling and minimal supersymmetry is the Klebanov- 
Strassler (KS) model (with an extension by Ouyang (4f| to incorporate fundamental 
matters). Although at the highest energies the gauge theory is best described in terms of 
bifundamental fields with effective degrees of freedom diverging, at the lowest energies the 
gauge theory resembles N = 1 SUSY QCD. Again in the limit when effective D3 brane 
charge is large, the gauge theory with large 'tHooft coupling has an equivalent description 
in terms of warped deformed cone. 

The thermodynamics of this non-conformal gauge theory is encoded in the dual geom- 
etry. A great deal of effort has been given in computing the black hole geometry dual to 
non-conformal thermal gauge theories. For example in [46, ^7], 48, 49] the cascading picture 
of the original KS model was extended to incorporate black-hole without any fundamental 
matter, while fundamental matter was accounted for in [50]. Most of the attempts are 
based on obtaining an effective lower dimensional action from KK reducing ten dimen- 
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sional supergravity action. Dimensional reduction of a generic ten dimensional action is a 
formidable challenge specially when there are non-trivial fluxes and scalar fields and it is 
highly non-trivial to obtain a consistent truncation. Furthermore, it is not clear how RG 
flow of the dual gauge theory can be obtained since the fields in the effective action are not 
the dilaton or the flux in the original ten dimensional action. 

On the other hand, partition function of the geometry and thus the thermodynamics 
of the dual gauge theory can be directly obtained by computing the on shell gravity action 



with appropriate boundary terms []51|. In a series of papers |52|- [pq| , we proposed the ten 
dimensional deformed resolved conifold black hole geometry as the dual to UV complete 
gauge theory which resembles large N thermal QCD. Working with the ten dimensional 
geometry, we avoid the difficulty of KK reduction while extracting the exact RG flow of 



the gauge theory. The thermal gauge theory studied in [52]-[pq] has a rich phase structure 
and as temperature is altered, we expect phase transitions. Our goal is to study the phase 
transitions at strong coupling, but to do so, we must first understand how geometries can 
describe different phases. In this work, we make progress in that direction and obtain 
the most general ten dimensional geometry (with or without a black hole) that arises 
from low energy limit of type IIB superstring theory. Directly identifying the gauge theory 
partition function with that of the ten dimensional geometry, we show how phase transition 
is realized. For any given temperature of the dual gauge theory, there are two geometries- 
extremal (without black hole) and non-extremal (with black hole) but the geometry with 
lower on shell action is preferred. At a critical temperature T c , both geometries are equally 
likely and we have a phase transition. 

Our ten dimesional solution is analytic and the corrections to the metric due to the 
black hole can be exactly written as a Taylor series in f^/r and g s M 2 /N where is the 
Schwarzchild horizon and M, N are number of five branes and effective number of three 
branes at some high energy. This series expansion allows us to write down exact expression 
for the on shell gravity action with or without the black hole and then obtain the critical 
horizon. However, in our analysis of the non-extremal geometry, we have considered con- 
stant axio-dilaton field. There are no D7 branes, no fundamental matter in the dual thermal 
gauge theory and no Baryochemical potential. On the other hand, it is straight forward to 
incorporate chemical potential in extremal ten dimensional geometry by considering gauge 
fluxes on holomorphically embedded D7 branes. Note that D7 brane embedding in black 
hole geometry is highly non-trivial in the presence of non-trivial three form fluxes 2 . What 
we have been able to obtain is the black hole geometry with only non-trivial five form F§ 
and three form G3 but constant dilaton and already in this simplified scenario, we find 
a first order phase transition. Observe that the warped deformed cone of KS model is 
dual to M = 1 SU(A4), M > 1 pure gauge theory which has no fundamental matter 
but it confines, exhibits a mass gap and has gluino condensates fjjjj. Thus the confine- 



2 In a particular scaling limit, the D7 embedding in a black hole geometry with considering gauge fluxes 
on the D7 branes along with UV completion was considered in J57|. The form of the localized sources was 
proposed but the exact value of the on shell action considering all the fluxes was not determined. Hence 
the exact value of critical temperature was not obtained, but the scaling of the critical temperature with 
the number of branes should be similar to the calculation done in this paper. 
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ment/deconfinement phase transition we obtain should mimic the phase transition between 
confined glue-ball and deconfined gluons of SU(N) QCD without flavor. In our upcoming 
work |B3] , we will incorporate the effect of running dilaton field and other localized sources 
on the non-extremal geometry which will also give the dual geometric description of UV 
complete thermal gauge theory. 

In section ^, we briefly discuss the supergravity equations and outline the procedure of 
obtaining thermodynamic state functions of the gauge theory. Summarizing the extremal 



geometry and its dual confined gauge theory in section 2.1 , we propose the non-extremal 
limit of ten dimensional geometry in section |2.2| . Using the supergravity solutions for 
extremal and non-extremal geometries, in section we study transitions between the 
geometries which essentially describes phase transitions in the dual gauge theory. When the 
boundary of the geometry r = 1Z scales with the 'tHooft coupling g s N, we find the critical 
temperature T c as a function of g s N. Finally in section ||, we discuss the structure of the 
gauge theory and propose a possible brane configuration that gives rise to the supergravity 
solutions. Since our ultimate goal is to learn about thermal phase transitions in nuclear 
matter, we describe the connections between the gauge theory and large N thermal QCD. 



2. Gravity Action and Gauge Theory 



We start with the type IIB supergravity action including local sources in ten dimensions 

mm: 
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in Einstein frame, where r = Co + ie~^ is the axio-dilaton with Go being the axion and 
(j) the dilaton field and F5 = (IC4 — \Ci A H3 + \B2 A F3 is the five-form flux sourced by 
the D3 or fractional D3 branes. Here G3 = F3 — tH^ with F3 the RR three form flux and 
H3 = (IB2 the NS-NS three form flux with B2 being the NS-NS two form while G2 is the 
RR two form. We also have G4 the four- form potential, G = \J — det gMN, M, N = 0, ..,9 
with gMN being the metric, i?( 2 ) the curvature two-form, and S\ oc is the action for localized 
sources in the system (i.e. D branes). 

The above action ( |2.lD is the most general supergravity action obtained from type IIB 
superstring action with fluxes and localized sources. Minimizing the action leads to the 
following Einstein equations 
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where [i,v = 0, ..,3; m, n = 4, ..,9 and we have assumed that the fluxes and axio-dilaton 
only depends on coordinates x m . 

The equation of motion for G3 can be expressed in terms of a seven- form A7 = *ioG3 — 
iCi A G3 in the following way, 

dA 7 - -Accfr A ReA 7 = (2.3) 
ImT 

where typically A7 quantifies the deviations from the imaginary self dual (ISD) behavior. 
Minimizing the action (2.1) also gives the Bianchi identity for the five- form flux |E| 

where p\ oc is the D3 charge density from the localized sources and T3 is the brane tension. 
In the presence of D3 branes, p\ oc is a delta function peaked at the location of the branes 
while D7 or fractional D5 branes may also contribute to D3 charge. 

To keep lorentz invariance along the space-time direction we assume the self-dual five- 
form has the form 

F 5 = (1 + *)da A dt A dx A dy A dz. (2.5) 

where a(x m ) is a scalar field, function of the internal coordinates x m . We take a general 
metric ansatz: 

ds 2 = _ e 2A+2B^2 + e 2A {(lx 2 + dy 2 + + e ~2A-2B g^m^n 

g mn dx m dx n = (a(r)dr + k(r)ds 2 M ^j (2.6) 

where B 7^ characterizes the existence of a black hole. Using this metric ansatz we can 
express A7 as 

A 7 = [e 4A+B * 6 G 3 - iaG 3 ] A dt A dx A dy A dz (2.7) 

The above choice of A7 leads us to three different classes of solutions from the G3 EOM 
Q2.3| ). These three classes can be tabulated in the following way: 

• If a = e AA+B in (|2.7j ) and A7 = dAj = then G3 must be imaginary self dual (ISD). 
When B = then this is the same as GKP solution [ 52 1 , and in this case r is not restricted 
and 

Rmn = !?7i d ~\ ~ 3V m d n B - h m Bd n B (2.8) 
2 I Imr \ z 2 

where R mn is the Ricci tensor for cone metric g mn - 

• If a ^ e 4A+B then we can take A7 7^ but keep cfA7 = and dr = 0. This means r is a 
constant 3 . 



3 Or r = d\-i i.e d of a (— l)-form. The functional form for the ( — l)-form is non-trivial, so this option 
is more cumbersome to use. 
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• If a 7^ e lA+B then we can again take A7 7^ but now dA? 7^ and dr 7^ such that (|2.5| ) 
is satisfied. This means both axion and the dilaton could run in this scenario. 



The first possibility is studied in the appendix [A|. We are more interested in the last 
two options because it is well established that AdS black hole solution corresponds to 
conformal SU(N) gauge theory, and non-AdS black hole solutions should reduce to AdS 
black hole solutions in certain limit. It is obvious that in the first case at G3 = it does 
not recover AdS black hole solution while the last two cases obviously do. 

Suppose now we find a solution to the above set of Einstein equations ( |2.2| ) along 
with the flux equations (2.3), (2.4). The resulting geometry with the metric ( |2.6| ) has the 
topology of y 5 x 5 5 and in particular for G3 = 0, Y 5 can be a five dimensional AdS space 
with or without black holes which has throat radius L while 5 5 can be a compact five 
sphere S 5 with radius L. For non- vanishing three form flux G3, Y 5 will be a non-AdS 
space and the size of 5 5 given by e~ 2A ~ 2B k(r) will diverge for large r. However for fixed 
radial coordinate r, the geometry has topology M4 x 5 5 where M4 is four dimensional 
Minkowski space and 5 5 has fixed radius. This way at fixed radial location (which we 
interpret as the boundary), we can obtain a four dimensional manifold M4 by integrating 
over the compact space S 5 . The gauge theory is defined on Minkowski space M4 while the 
manifold Y 5 x 5 5 is the dual geometry for the gauge theory. 

According to AdS/CFT correspondence, the partition function of string theory on 
AdS$ x S 5 should coincide with the partition function of M = 4 super- Yang- Mills theory. 
For generalized gauge/gravity correspondence, we can identify the gauge theory partition 
function with that of ten dimensional dual geometry, 

7 — p~ F / T —7 ■. ~ p^CTavity 

-■^gauge — c — -'-'gravity — e 6 ' 

'-'gravity = "^total + &GH + 'S'counter (2-9) 

where F,T are free energy and temperature of the gauge theory, Stotai is given by (^Q|), 
Sqh is the Gibbons-Hawking boundary term [51] and S counte r is the counter term necessary 
to renormalize the action. The gravity action is evaluated on Y 5 x 5 5 and then wick rotated 
t = it to obtain the Euclidean on shell value. The ~ is because we have ignored all the a' 
corrections and loop corrections. 

Using ( |2.9D , we can in principle obtain all the thermodynamic quantities for the four 
dimensional gauge theory by considering the ten dimensional action 5*g^ vity on shell. For 
instance, free energy F and internal energy E of the gauge theory are given by 

F = -T log (£ gaugc ) = -T log (^gravity) = T S^ vity 
n a oren 

E = T^(lo g Z^ gc ) = ^L (2 . 10 ) 

where T is the temperature. Knowing the free energy one gets the pressure p and entropy 
s 

fOF' 



where V3 = J d?x is the volume of three dimensional flat space. 

Our primary concern for this paper is to obtain the ten dimensional geometries that 
can arise from the action ([O]). However, the same action can give rise to more than one 
manifolds. In fact, just like the case for AdS space discussed by Hawking and Page [63] 



and elaborated by Witten Q , there are two manifolds X 1 = S 1 x M 1 and X 2 = S 1 x M 2 
that minimizes the action (^T]). The manifold which has lower value for the on shell action 
for a given temperature of the dual gauge theory will be preferred. Since X 1 and X 2 are 
distinct geometries, the thermodynamics of the gauge theory will be different at different 
temperatures, depending on which geometry is preferred. This means X 1 and X 2 will 
correspond to different phases of the gauge theory and we will now analyze the manifolds 
in some details bellow. 

2.1 Extremal geometry and confinement 

The metric of the extremal geometry X 1 without any black hole, i.e. B = 0, is given by 
Ejj] (with Minkowski signature) 



ds 



2 _ J2A 



dt 2 + dx 2 + dy 2 + dz 2 



+ e- 2A g mn dx m dx n (2.12) 



where e iA = h{p) and g mn is the metric of the deformed cone 

g mn dx m dx n = l -A^K{p) {dp 2 + ( 5 5 ) 2 ) + cosh 2 (0 [(g 3 ) 2 + {g'f 



i 2 (f ) [(a 1 ) 2 + {a 2 ? 

where A is a constant, g l ,i = 1, .., 5 are one forms given by 



+sinh 2 (J) [(5 1 ) 2 + (9 2 ) 2 ] \ (2-13) 



1 e 1 - e 3 2 e 2 - e 4 

9 = ^7T< 9 = ^JT 

e 1 = — sin#i d^i, e 2 = 

e 3 = cost/> sin#2 ^2 — simp d9 2 , 

e 4 = smi/j sm02 dfa + cos^ d6 2 , 

e 5 = d^ + cos^i d0i + cos6>2 #2 (2.14) 



and 



The three form fluxes G3 on the deformed cone is Imaginary Self Dual (ISD) to preserve 
the M = 1 supersymmetry while the five form fluxes is self dual 

* 6 G 3 = i G 3 , F 5 = (l + *)dh~ l Adt Adx Ady Adz = dC A + B 2 A F 3 . (2.16) 
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Note that the metric ( |2.13 ) is Ricci flat, that is the Ricci tensor for the metric g n 



denoted by R mn = (but Ricci tensor for the metric g mn denoted by R mn / 0), while 
minimizing the action ( |2.1| ) with non-zero dr and localized sources will give rise to an 
internal metric which is not Ricci flat in general. In particular, the Ricci tensor for the 
warped metric is 

Rfw = ^;^77yV 2 iog/i(p) 

Rmn = Rmn ~ V 2 log/l(/9) - -d m logh(p)d n logh(p) (2.17) 



where 



V 2 = g mn d m d n + d m g mn d n + \g mn g pq d n g pq d m (2.18) 



Now using ( p. 17 ) in ( |2.2j ) with the metric given by ( |2.12 ), one readily gets that R 



dmjd m T + local sources. However, F-theory |6(], |6?J dictates that the running of the axio- 
dilation field is dr ~ 0{g s Nf) where Nj is the number of seven branes and we will only 
consider terms which are of 0(g s Nf) ignoring all higher order terms. If we also ignore the 



localized sources, then indeed we obtain R mn = f57], p2;] . 

Even though g mn remains the same the warp factor is different in different regions. In 
the IR region we assume the D7 branes are far away and thus the axion-dilaton is constant. 
The warp factor, B 2 and F3 in this region are given in [Q] and calculated to the first order. 
The warp factor can be written in the following form 

h(p) = (Hp* (2.19) 

where the coefficient Cj can be treated as constants. When doing so it is assumed that 

*dh~ 1 A (ft A dx 1 A dx 2 A fix 3 = B 2 A F 3 (2.20) 

This is correct if the number of D3 branes is multiple of that of the D5 branes. In the rest 
of the paper we will restrict to this special case 4 . 

The warped deformed conifold with h(p) given by ( 2.19| ) for small p as proposed in 



1 44] in fact removes the IR singularity of the warped conifold |35| and gives rise to linear 
confinement in the dual gauge theory. Unlike the regular cone, the deformed cone has a 
blown up S 3 at the tip of the cone and this finite size of the S 3 removes the IR divergence 



of the fluxes [44]. Observe that near p = 0, the metric in ( |2.32j ) reduces to that of S 3 

1 



l dx m dx n ~ ^ 4 / 3 (2/3) 1 / 3 



(2.22) 



4 In general 



A dt A dx 1 A dx 2 A dx 3 = B 2 A F 3 + Ng 1 A g 2 A g 3 A g 4 A g 5 , (2.21) 



where is the D3 brane charge at certain energy scale ro where B2AF3 = 0. If is written as N = IM + k 
with < k < M, then IM can be absorbed into B2 A F3 leaving only the k term. This is the cascading 
shown in the gravity side. 
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which implies that the radius of S 3 at the tip of the cone is of 0(A 2/3 ). The finite size of 
S" 3 at p = gives finite value for the three form flux strength I-F3I 2 , l-f^l 2 at the tip of the 
cone. This way, the IR singularity of the fluxes are removed, which consequently removes 
the IR singularity of the warp factor. 

On the other hand, this constant A modifies the embedding equation for the cone 
and is related to the expectation values of gauge invariant operators 5 in the dual gauge 
theory H]. Hence a non-zero A results in non-zero expectation value and corresponds 
to spontaneous breaking of Z<im group down to Z2 where M is the number of fractional 
three branes placed at the tip of the cone. Note that in the absence of D7 branes, this 
spontaneous broken symmetry is not identical to the chiral symmetry of QCD as there are 
no fundamental flavor. But even in the presence of seven branes, the deformed cone metric 
Q2.12 ) with warp factor ( p. 19 ) will be valid by considering the D7 branes as probes. The 
probe seven branes will give rise to fundamental matter and depending on their embedding, 
can lead to breaking of flavor chiral symmetry just like in QCD ]68j ] . 

Coming back to the metric ( |2.13| ), observe that with a change of coordinates 



A 2 e p 



(2.23) 



for large p, the metric becomes 

g mn dx m dx n ~ dr 2 + r 2 I ~ 



i=i / 



(2.24) 



which is the metric of regular cone with base T > . Thus only for small radial coordinate 
p, the internal metric is a deformed cone while at large p, we really have a regular cone 
with topology of R x T 1 ' . 

However for large p, we can no longer ignore the running of the r field as we will 
be near the seven branes. As we have a regular cone for large p, we can use Ouyang's 
holomorphic embedding of seven branes to determine the running of the r field along 
with the modified flux G3 (which is again ISD) and -F5. Then for p > 0, using change of 
coordinates ( p. 23 ), we obtain the warp factor 



+ 



h(r,e 1 ,e 2 ) -- 

81g 2 M 2 N f 
32^ 



a' 2 I 

7T «o + 



81g s M 2 



log - 



, , 3ff s iV> , 
1 + ^ (1 ° g 



log (i) log ( si 4 si 4) 



1 

+ 2 



(2.25) 



where Nj is the number of D7 branes which is present in the action (2.1) and source the 
r field while r m ; n is the minimum radial distance reached by the seven brane. The exact 
brane configuration that could give rise to such a warp factor in the dual geometry will be 
discussed in section |||. Also ao, r\ are constants which will be determined by matching the 



5 At the bottom of the cascade when all D3 branes disappear and only fractional branes remain, we have 
pure glue SU(M) with gluino condensates (AA) Q. The dual geometry partially captures the feature of 
the pure glue theory and we expect non-zero A corresponds to non-zero (AA). 



-9- 



above solution for large p with the solution ( |2,19| ) valid at small p. In fact ao is proportional 
to the D3 brane charge at r = 77. To see this we notice that on the gravity action, we do 
not have any branes - only fluxes. However, due to the duality, the units of the fluxes on 
the gravity side should equal to the number of branes on the gauge theory side. This way 
the units of five form fluxes can be identified as the effective charge: 

^efr = ^r/ h = -^r^ /r^WA/A/A/ (2.26) 
2k( T 3 J T i,i 216K| r 3 J dr 

where T3 = ^3 = (2ir)~ 3 a'~ 2 is the three brane tension and the Gaussian surface is warped 
T 1 ' 1 with radius y/ r 2 \/h that encloses Minkowski space. Using ( |2.25 ) in ( 2.16 ), we readily 
get from ( [2.26 ) that N e g depends on r. In particular at r = n 



N cS (r = n) = §^ + 0(g s M 2 ,g 2 M 2 Nf) 

V 5 = J d^d^idfacteM sin6»isin6» 2 (2.27) 



Using the exact same argument, we get N e s(p = 0) = as dh/dp = at p = 
This implies that there are no effective D3 branes, they have cascaded away at the IR and 
we are left with only fractional D3 branes. 

Note that the holomorphic D7 brane embedding of Ouyang jl5| only gives U(Nt) 
flavor symmetry group. However by considering Nf number of D7 and D7 branes, the 



chiral Ul{Nj) x Ur(NA symmetry can be realized |68|, ^] but then the warp factor ( gjgp 



will also be modified. In addition, the warp factor ( |2.25| ) only makes sense in the large 
p region which means r m ; n is large. Thus the fundamental matter fields resulting from 
the holomorphic embedding have large mass and we do not expect chiral symmetry to 
hold for massive flavors. As supersymmetry is not broken at zero temperature for the 
Ouyang embedding, the warp factor ( p. 25 ) corresponds to a gauge theory with massive 



fundamental matter that is super symmetric at zero temperature. Recently in []69[| , the 
dual to non-supersymmetric confining gauge theory was proposed and the metric ( |2.12 ) is 



consistent with this proposal with the understanding that g mn in ( 2.32; ) be replaced with the 



perturbed metric g} nn which is also Ricci flat. We will analyze possible brane configurations 
that may give rise to SUSY and non-SUSY gauge theories with fundamental representation 
matter in section ^. . 

The solution ( [2.25 ) leads to divergence of N e g for very large r i.e very large energies. 



But this is expected as the dual theory is a cascading gauge theory with effective color 



diverging in the UV [44, 59, 70]. However our goal is to study thermal QCD which becomes 
conformal and free in asymptotically large energies. Thus the dual geometry with warp 
factor ( p. 25 ) cannot be relevant for a QCD like theory and must be modified for r ~ A — >• 00. 



In [|54], [58|] a proposal was made to modify the UV dynamics of the cascading gauge theory by 
introducing anti five branes separated from M D5 and D3 branes along with embedding 
D7 — D7 branes to account for fundamental matter and chiral symmetry breaking. The 
additional branes modify G3 and F5 while the form of the r can be obtained from F-theory 
1 66 , 67]. Note that the D5 — D5, D7 — D7 system behaves as dipoles and thus away 
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from branes, their contribution to flux and axio-dilaton is not divergent and behaves as 
l/r l ,i > 0. One can consider brane embeddings such that the fluxes and axio-dilaton fields 
in the dual geometry is well defined everywhere in the bulk geometry while the form of the 
metric ( p. 12 ) remains unchanged. However, the internal metric does not remain Ricci flat 
anymore and the metric (2.24) gets corrections of O(ggNj). 

The fluxes proposed in [54| essentially give B 2 — > for r > ro 6 while the axio 
dilaton field behaves as r ~ bi/r l , bi ~ 0(g s Nf) for large r. Thus for asymptotically large 
distances, we end up with constant axio-dilaton field and a vanishing NS-NS two form B 2 
and the geometry behaves as an AdS§ x T ' . This means that addition of anti five branes 
have slowed down the cascade and the theory has no Landau poles and the associated UV 
divergences. The warp factor becomes 



for large r where O = <fo, ip). Using (|2.26| ) one readily gets that N e s(r = 00) ~ ao 
and hence does not diverge. As N e Q reaches a fixed value with the geometry becoming 
AdS§ x T 1 ' , the gauge theory reaches a conformal fixed point. However, for supergravity 
approximation to be valid, we need A fj(oo) ~ ao 3> 1 which means the 'tHooft coupling 
g s N c fi(co) is still very large . Although we have a conformal field theory at large scale 
A = r — > 00, the gauge group has large number of effective colors with large 'tHooft 
coupling and thus the gauge theory is not asymptotically free. However, the Yang Mills 
coupling of our gauge theory (51,02) behaves the following way 



1 


1 




+ — 




92 


1 


1 




~7 2 



B 2 (2.29) 

s 2 

At asymptotically large energies, that is r — > 00, B2 — > and we have 01 — 02 = gYM with 

2 1 . . 

lim = — (2.30) 

A ^°° 9ym 9s 

Now we can take the limit g s — > by keeping g s N e g fixed, which means gyu — > at 
asymptotically large energies, just like in QCD. However, since g s N e g is held fixed at large 
values, we have large 'tHooft coupling while for QCD, 'tHooft coupling goes to zero. 

In summary, the extremal geometry takes the simple form given by Q2.12D while the 



warp factor is given by (|2.19|,|2.25|) and ( f^g) for various range of the radial coordinate p. 



The fluxes for various regions of the geometry can be found in [p4|]. 

Now to obtain thermodynamic state functions such as free energy, entropy and pressure 
of the gauge theory that is dual to the geometry X 1 = Y 5 x5 5 , we must obtain the on 



6 ro ~ Ao is related to the higgs mass in the dual gauge theory. For details of the higgising mechanism, 
please consult |B8| 
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shell supergravity action on the manifold X 1 . However note that, since there is no black 
hole horizon in the geometry X 1 , the Euclidean renormalized on shell action Suavity = 



is independent of horizon with F independent of T. Then using (2.10) one obtains that 



the free energy is independent of temperature and ( |2.11 ) readily gives 

dF 

* = -^ = (2.31) 

The above result for entropy is consistent with the confined phase of large N gauge 
theory. For confined phase, Zn symmetry is preserved and physical quantities do not 
depend on the temporal radius at 0(N 2 ), due to large iV volume independence f7lf , |fr2| . 
This means free energy is independent of temperature and entropy is zero at 0(N 2 ). Thus 
the gauge theory dual to the geometry X 1 is in the confined phase. We will later see that 
for X 2 , entropy is nonzero- indicating we have a deconfined phase and degrees of freedom 
have been released at high temperature. 

2.2 Non extremal geometry 

The non-extremal geometry X 2 is a regular black hole with B ^ in (^1]). The internal 
manifold is typically a resolved-deformed conifold that mimics the symmetries (or breaking 
of symmetries) of the dual gauge theory. For simplicity, we will consider constant axio- 
dilaton field, that is there are no seven branes in the dual gauge theory, i.e Nf = 0. We 
make the following ansatz for the internal metric 

g mn dx m dx n = a{r)dr 2 + k{r)g$ dx p dx q 



A 



~g^dx p dx q 



+ ^(E^) 2 )J + 9' pq dx p dx q (2.32) 

i=l 



where g l are defined in ( 2.14| ), k(r) = r 2 e 2B b(r) and x p ,x q represent angular coordinates. 



Note that when there are no black holes i.e. horizon = and no three form flux i.e. 
M = 0, we have g' pq = and a(r) = b(r) = e 2B = 1. On the other hand, as we will see 
explicitly, fluxes enter the Einstein equations with terms of 0(g s M 2 /N), where 



54K 2 n T 3 N eS (r = rA 

N = a ~ 10 t/, 9 V 2.33 

V 5 a' 2 



is proportional effective number of D3 branes at certain energy scale 77. These scalings 
with horizon and flux strength indicate that at the lowest order, we have 

9' pq = 0(r h9s M 2 /N) 

a(r) = 1 + 01 (r) = 1 + 0(r h g s M 2 /N), b(r) = 1 + 6 x (r) = 1 + 0{r h g s M 2 /N) 

(2.34) 

Note that with g' pq 7^ the two S 2 are squashed and topologically the base of the cone is 
no longer S 2 x S 5 . 

Before solving the equations of motions, we would like to discuss some subtleties: 
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• We are studying the gravity dual of the gauge theory and these gravity solutions are 
generated by fluxes only. The gauge theory lives on N D3 branes and M D5 branes placed 
at the tip of a regular conifold, the back reactions of the branes will give us a warped 
geometry. There are fluxes in these gravity solutions whose integral J F is non- vanishing, 
however that does not mean there is a S function or brane source in such solutions. Actually 
J F can be non-zero without any source if the volume of integrated space is non- vanishing 
anywhere. This is what happens in AdS§ x S 5 case, J s5 F§ ^ where S 5 is a 5D sphere with 
constant radius. Particularly in Klebanov-Tseytlin solution, the radius of the warped T 1 ' 1 
becomes zero at non-zero r and becomes ill defined at very small r, leading to singularities. 



This singularity was resolved by Klebanov and Strassler [44| who pointed out that at 



low energies, the gauge theory superpotential receives quantum corrections [T3] and the 
corresponding dual geometry is no longer the warped regular cone but it must be warped 
deformed cone whose S 3 is non- vanishing at r = 0. Thus even though J TM F5 7^ and 
J$3 F3 7^ there is no 5 function in either flux equations or Einstein equations. In section 
(Q) , we will outline the brane configuration that may give rise to the dual gravity solutions 
that we present here. 

• The Schwarzschild black hole is a vacuum solution of Einstein equations with the stress 
energy tensor T mn = 0. The horizon radius is an independent parameter and not deter- 
mined by the Einstein equations. By considering the weak field limit, one can obtain the 
geometric mass of the black hole from the horizon radius which is a free parameter. Simi- 
larly by placing D3 branes at the tip of regular cone, we can generate AdS Schwarzschild 
black holes and again the black hole mass is a free parameter independent of the number 
of D3 branes. In our solution discussed below we will find a charged black hole and we 
interpret the charge and mass arising from the D branes placed at the tip of the regular 
cone. We will see that the black hole is not a Schwarzschild black hole and the horizon 
is no longer independent of the number of D branes. However, the solution can be built 
from the Schwarzschild solution and thus the vacuum solution plays a crucial role in our 
analysis. 

Now lets look at the equations of motions. We first consider the external components 
of the Einstein equations in (^Z 
becomes 

G3 • G3 



Using the form of -F5 (|2.5| ), the first equation in 



R 



48 Imf 



+ 



+ K10 



T, 



loc 



fj,l> 



loc 



(2.35) 



On the other hand, the Ricci tensor in the Minkowski direction takes the following simple 
form 

1 r 



R 



2 L 



d m (g mn d n g^ l/ ) + g mn T^f M d m g^ u 



g mn g v ^d m9ll ^d n g v 



where v 1 , fi' 



0, .., 3 and T^f M is the Christoffel symbol. Now using the ansatz 



(2.36) 
\j for the 



metric, (2.36) can be written as 



Rtt = e^ A+B ^ [V 2 (A + B)- Zg mn d n Bd m {A + B)\ 



(2.37) 
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where the Laplacian is defined in (|2.18|) , 

The set of equations can be simplified by taking the trace of the first equation in (| 
and using ( 2,37] ). Doing this we get 

^ {AA + B) _ 3g mn d n Bd m (4A + B) = e -2A-2B^^_ + ^ 10A _ AB 

12Imr 

+ ^M e -2A-2B ( jnn _ T ploc (2.38) 

On the other hand using ( |2~35D in ( |2~3?1 ), one gets 

Rl-R x x = (2.39) 
which in turn would immediately imply 

V 2 S - S^dmBdnB = (2.40) 
Now let's look at the Bianchi identity. Using fl2.5|) in (|2.4[) gives 



V* a _ 3e- 2A ~ 2B d m Bd m a = 6 ^-b HG 3 -G 3 + 2e _ 6A _ 3B g^A+B gm a (2 41) 
Subtracting it from ( [2.38 ) one gets the following 

v2( e 4A+B _ a) = ^l\ iG3 _ * 6 G 3 |2 + e -GA-3B ld{e 4A+B _ q) |2 

24Imr 

+3e- 2A - 2B a mJ Ba m (e 4A+B - a) + local source (2.42) 

Next we consider the internal components of Einstein equations. Ricci tensor for the 
x m , m = 4, .., 9 directions take the following form 

Rmn = Rmn + fan V 2 (A + B) - 3g mn g Xk d x Bd k (A + B) 

+ 3V m d n B + d m B8 n B - 8d m Ad n A - 2d (m Ad n) B (2.43) 

where V m is the covariant derivative given by 

y m V c = 8 m V c - f b mc V b (2.44) 

for any vector Here R m n is the Ricci tensor and V^ nc is the Christoffel symbol for the 
metric g m n- 

Using (|2~38|) and ( g!4C| ) we find 

~ _ G3 " ^3 , G ma i, • G^ 6 d m Td n T 

-t<"mn — fan n . T ~ T ~ t 



24Imr 4Imr 2 | Imr \ 2 

_ l e SA-2B dmadna + 8dmAdnA 

-3V m 8 n B - 8 m Bd n B + 2d [m Ad n) B (2.45) 

To solve for these equations we must first find the three-form flux. For simplicity we 
will ignore the running of the r field, that is dr = 0. We further take Co = and = g s . 
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With constant r and considering fluxes on the deformed cone, G 3 = F 3 — tH 3 must be 



closed everywhere since we are ignoring all localized sources i.e. S\ oc = in (2.1). Then 
the three form flux equation becomes 

d(e 4A+B * 6 G 3 - iaG 3 ) = (2.46) 

Written in terms of F 3 and H 3 , the above equation becomes 

dCi = 0, Ci = e AA+B * 6 F 3 - ae~*H 3 

dC 2 = 0, C 2 = e AA+B * 6 e^H 3 + aF 3 (2.47) 
where C\ and C 2 are related by 

/ 8A+2B \ p 4A+B 

* 6 C 2 = a *6*3 + Ci. (2.48) 

\ a / a 

Our goal is to solve the system of all Einstein equation and field equations exactly. 
To do this we take the simplest case where we assume all the unknown functions A, B, a, 
b and a only depends on radial direction r. In the non-AdS extremal limit and the AdS 
non-extremal limit, all these functions are known, so we can use the known forms as a 
starting point. We use the following ansatz for the warp factors 

iA i / \ Na' 2 ( 81o s M 2 , (t\ , . , 

e = h ^ = — { 1 + Hlv- log (n) + Mr > 

e 2B = g(r) = 1 - ^ + G(r) (2.49) 

where 77 is some scale determining boundary condition of h(r) and fh is the Schwarzschild 
horizon i.e. the black hole horizon in the limit M = 0. We will later discuss the relation 
between the actual horizon rh and fh- But for now note that for F 3 to recover the Klebanov- 
Tseytlin solution in the limit B = 0, the corrections must be proportional to the black hole 
horizon r^. On the other hand, these corrections must vanish as r 1— > 00. However F 3 must 
not depend on r since D5 brane is orthogonal to the radial direction. Thus we argue that 
F 3 remains the same 

Ma' 



-g b A oj 2 (2.50) 



2 

where uo 2 = g 1 A g 2 + g 3 A g 4 . On the other hand, using this form in Q2.47| ) gives 

aH3 = (^^l ia,-*,) (2.5D 



where we have used the ansatz ( 2.34| ),( p.49| ) and ignored terms of order O (rf l g s M 3 /N 



2) 



which is a good approximation in the limit 3> M. If we assume C\ = c\(r)dr A uj 2 , this 
automatically gives dH 3 = 0. Now from the closure of C 2 and using ( p.48| ), we find 
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so we get 



e-*H, = a ^ (2.53) 

We further compute G3 • G 3 

ri s=i j 6 ft A *6G3 

G 3 • G 3 dx u = — — 

= ^W-(l + (a - q^))^ 5/1 ^) sin 0i sin 02/^56^ 



\/<76 = 



e~ 6A ' B r 5 Vabb 2 (1 + 0{r h g s M 2 /N)) sin0i sin0 2 



108 

so 



01 M 2 n l2 e %A 

G 3 -G 3 = 81M r6 a b3 6 (1 + (a - a(r,)) 2 e~ 8A - 2B ) (2.55) 

where we have ignored terms of O (r^M 3 /N 2 ) that arise due to the deformation of the 
base g'pq- 

On the other hand, Gm^G^ is 

^mab^n — r mab r n ' e ^mab^n 

= ^(*F 3 ) ■ (*F 3 ) - (*F) mab (*F)f + e- 2 *H mab H« b 
17M 2 n' 2 p &A 

= 9mn r%3mf ~ (*F) mab (*F)^ + e- 2 *H mab Hf (2.56) 

While T mn = -(*i ? ) ma fe(*-? 1 )^ + e~ 2( f'H rnab H^ b only has following non-zero components 

4 2 a' 2 e- 2B a ( 8A _ 2B _ 2 _ x 72 

162M 2 a' 2 e 6A ( RA 9R/ , \ 
9MVV^«/ -%a-2B, f ^2 A 6sin 2 0, 

243M 2 a ,2 e 6A sin 2 0; / _ 8A _2 B/ , Nx2 n \ 

81M 2 a /2 e 6A / 84 2B/ , ,., \ 
= ^ ^53 (e- 8A - 2B (« " «(^)) 2 " l) (2-57) 

where again we have ignored terms of O (r^^M 3 /^ 2 ) that appear due to the deformation 
of metric of base g 1 . 



9M 2 < 

T — 



= grr- 

T <t>i<S>i = 
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Using the above forms of the fluxes, the only non-trivial internal Einstein equations 7 



( |2^5D are 



Rr 



297Q 1 



-SA-2B g a g a + 8 Q rA Q rA 



R 



-3S7 r d r B 
27Q 



- d r Bd r B + 4<9 r Ad r B 
+ 3f r ^d r B 



Rj 



Rr! 



~9<t>l 



27Q 



Rfi 



27(15 sin 2 0* - 2 cos 2 6 t )Q 
1 8(2 cos 2 6>; + 3 sin 2 (9;) 



+ 3r!l, ^ct-B 



(2.58) 



where Q = M 2 a /2 e 6A (e" 8A_2B (a - a(r fe )) 2 - l)/(r 6 6 3 Imr) and j = if), fa. 

We will now solve the set of Einstein equations ( 2.38[ ),( pT4C ) and (2.55) along with 



the trivial Einstein equations and the Bianchi identity ( 2.42 ). Note that with the ansatz 
( 2.32 ) for the metric, the three form fluxes H3,F^ which satisfy the flux equation (|2.46| ) 
are already determined and given by (§~J0|) and (|2~53|) . Thus, once equations (|2.38|) , (^40|) , 



(2.42) and (2.58) are solved, we have a complete set of solution as long as the trivial 
Einstein equations are also solved. To proceed with the analysis, we make the following 
ansatz for the five form flux strength 



t(r)e 



AA 



t{r) = 1 + h(r) = 1 + 0(r h g s M 2 /N) 



(2.59) 



which is again consistent with the extremal limit (i.e. = 0) and the AdS limit M = 0. 
The 0(rhg s M 2 /N) term is exactly consistent with the ansatz for the warp factors (|2~49| ) 
and the internal metric ( |2.32 ), (2.34). 

Also observe that the above ansatz (2.59) along with the warp factors (2.4S) reduce 
to the well known Klebanov-Witten black hole solutions in the limit M = 0. In fact, 
for M = 0, the horizon is given by with A\(r) = G(r) = in(2.49) and we have a 
Schwarzschild black hole in AdS*, x S 5 . We will refer to r~h as the Schwarzschild horizon 
while r/j denotes the black hole horizon in the presence of three form flux and other sources. 
With no black hole but M ^ 0, our ansatz reduces to the extremal limit of Klebanov- 
Tseytlin model with 7% = = B = A\ = G = 0. With these limiting behavior and the 
scaling of the warp factors, it is reasonable that the functions A\{r),G{r) appearing in 



(2.49) should have the following scaling 

A 1 = 0(r h g s M 2 /N) 
G = 0(r h g s M 2 /N) 



(2.60) 



With the above ansatz for the warp factors and ft, the system of equations can be 
further simplified. Using the expansions ( |2.49| ), ( 2.59|) and ( |2.60 ), we find for r ^ , 



7 By trivial Einstein equations we mean Ricci tensor, R m n = equations. 



-17- 



expanding up to linear order in g s M 2 /N 



~ i 



e - 8A - 2B (a - a(r h )) 2 = 1 - -£ + 0(r h g s M 2 /N) (2.61) 

Now the Taylor series in f^/r away from r = f/j 9 , keeping only up to linear order terms in 
0(g s M 2 /N) gives 



~4 

g pq Q ~ ^|0( 5s M 2 /iV) 

' J, 



~ E ^fe°(^M 2 /iV) (2.62) 
fe=i 



But in the limit N S> M, we have 0(g s M 2 /N) <C 1 and for f^/r < 1, we can ignore g mn Q 



all together. This drastically reduces equations in (|2.58|) to 

R rr = -^e- 8A - 2B d r ad r a + 8d r AdrA 

-3V r d r B - d r Bd r B + 4d r Ad r B 
Rij = 3T^jd r B 



(2.63) 



In fact, the above equation simplifies further, with the following observation: in the limit 
M = 0, ( |2.63| ) is trivially satisfied for our ansatz of the warp factor ( p.49| ) along with 
the scalings ( p . 60 ) , ( 2T5S| ) and ( |2.34| ). For M / 0, the right side of the equation has 



terms of order 0[ff /r 4k )0(g s M 2 /N), if we ignore the deformation of the base (that is 
a = b = 1,G = and g' pq = 0) and consider t\ = A\ = 0. On the other hand, the left 
side Rmn can have terms of 0(f 4+l /r 4+l )0(g s M 2 /N) if the base of the cone along with 
t\,A\,G have terms of same order. This means, to solve ( [2.58 ) everywhere from horizon 



to boundary along with the trivial internal Einstein equations, we expect the following 
scalings : 



oo /- \ 4+1 



1=0 V 7 V 



-~f:(*f"(*£'i 

1=0 v 7 v 



where T\ = Ai(r),ti(r),ai(r),bi(r) or G(r). Away from the Schwarzschild horizon i.e 
r > f h such that 0(g s M 2 /N) ~ 0(rfJr A ) < 1, we can ignore 0{g s M 2 /N)(D(rfJr 4 ) and 
then ( p. 58 ) is trivially solved. In fact, ignoring Ax(r),ti(r), ai(r), b\(r) and G(r) away from 



8 The Taylor series in g B M 2 /N is well defined Vr except r = fh and thus the expansion is valid away 
from but arbitrarily close to r = fh- 

9 The Taylor series for fh/r is well defined even at r = fh, for M =fc 0. 
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the horizon gives that all the equations ( |2 . 38| ) , (|2 . 40[) , ( |2~l38| ) along with the Bianchi identity 



(2.42) are satisfied with the following solution 

-4A h f ^ Na ' 2 (a j_ 81 9sM\ r 
e = h (r) = — ^ 1 + log- 

e 2S = 50 (r) = 1 - ^ ( 2 - 65 ) 

Now what is the relation between the real horizon and the Schwarzschild horizon 
r/i? To answer this, note that since A?" 3> M, we do not expect G to be very large even at 
the horizon. This is also consistent with the scaling in ( [2.64 ). This means, the true horizon 



rh is close to the Schwarzschild horizon f/j. More precisely we expect 



10 



r h = f h (1 + 0(g s M 2 /N)) (2.66) 



Thus r > fh also means r > and hence the warp factors ( |2.65 ) with ignoring 
•FiiQ'pq works exactly away from the horizon. At the horizon, the Ricci tensor equations 



( 2.63 ) are not valid and we have to incorporate the deformation of the base metric g' pq 
in solving the system of equations. In fact exactly at the horizon we can see that g mn Q 
is no longer ignorable up to linear order since fh/f ~ 1. If g' pq = and we consider the 
near horizon region, there are more than two independent equations in ( |2.58| )- resulting 
in more equations than the number of independent functions i.e Ai,G,ti,ai,b±. Thus we 
must consider the deformations g' g near the horizon to solve all the Einstein equations in 



(2.58) and the flux equations exactly. The form of the ansatz (2.64) is precisely chosen to 
solve all the equations in ( |2.58| ) along with the trivial equations. 

As long as we are away from the horizon, the details of the metric g' pq and functions 
Ai(r),ti(r),ai(r),bi(r),G(r) will not effect our calculations since all these functions are 
second or higher order in 7 ~ 0(g s M 2 /N) ~ 0(r£/r 4 ). Note that the flux equation ( ggjj ) 



is solved exactly in the entire region < r < 00 up to linear order in 0(rhg s M /N) 
and this equation does not involve g' pq . On the other hand, using (Hp in djHD, one 
would obtain that #3 will get contribution from g' . But what enters into the system of 
equations are not just H3 but the tensor T mn and gmnGz ■ G3, which up to linear order and 
away from the horizon give rise to the simplified equation ( 2.63| )- which does not involve 

9 pq- 

In summary, we have been able to obtain the exact form of flux and metric in the 
presence of a black hole. The three form flux is given by ( 2.5Cp and ( 2.53j ), while the scalar 



function that gives rise to the five form flux is given in ( 2.59| ). The modification t\ to a 
from the extremal limit, the modification of the warp factor A\,G and deviation of the 
internal metric is described by functions a\,b\ and g' pq - all of which take the form ( 2.64| ). 



Thus if rh > fh, we have exact expressions for the metric and the fluxes in the entire region 
r h < r < 00, up to linear order in 0(g s M 2 /N). Most importantly, using these solutions 



10 We expect D5 branes to increase the mass of the black hole and thus rh > fh- This is also consistent 
with our numerical analysis in J57J. 



- 19 - 



and the solution for the warp factor ( |2.49D , one can exactly compute the three form flux 
strength up to linear order in 0(g s M 2 /N)- a result that allows one to compute the on shell 
action which will be shown in section |3,2.1| , 

3. Phase Transitions 

Phase transitions of four dimensional SU(N) gauge theory at finite temperature can be 
realized by spontaneous breaking of the center symmetry Zjy along the Euclidian time 
circle. The order parameter associated with this symmetry is the temporal Polyakov loop 
operator 



1 



W = -TV P exp ( i I dr\ a A a ) (3.1) 



o 



where /3 is the periodicity of Euclidian time, P denotes time ordering and A a are the Gell- 
Mann matrices for SU(N). In the confined phase, Zfq symmetry is preserved and (Wo) is 
zero. In the deconfined phase, Zn symmetry is spontaneously broken and (Wo) ^ 0. 

In [54], we computed (Wo) in both extremal and non-extremal geometries following 



the procedure of p4|]-|77||. We found that for the extremal geometry, (Wo) = while for 
the black hole geometry, (Wo) ^ 0. Although the calculation was done in UV complete 
geometry, it is straight forward to demonstrate that for a Klebanov-Strassler geometry with 
a deformed cone and constant warp factor h(p) ~ constant ^ near p ~ 0, the Nambu- 
Goto action Sng ~ P-^qq ~ i- e - ^ linearly increases with inter quark separation T>, 
for large T>. Since (WoWj) ~ exp(-/3F Q g) and the free energy of a single quark is obtained 
in the limit Fq = limx>-->oo we get (Wo) ~ exp(— [3Fq) = 0. For the non-extremal 

geometry, even with non-UV complete scenarios, one can follow the procedure in and 
obtain that (W ) / 0. 

Thus the extremal geometry X 1 corresponds to zero order parameter while black hole 
geometry X 2 describes non-zero order parameter. Furthermore, as already discussed in 
section 2.1, using ( |2.11 ) one directly obtains that X 1 corresponds to zero entropy while 



X 2 has black holes and non-zero entropy. Thus X 1 is the gravity dual to the confined 
phase of the gauge theory while X 2 describes the deconfined phase. Note that although 
X 1 corresponds to zero entropy, it can describe the gauge theory at non-zero temperature 
and we can obtain the free energy of the dual gauge theory by using ( 2.10[ ) with the on 



shell action for X 1 . This is consistent with the understanding that confined phase can exist 
not only at zero temperature, but up to some non-zero critical temperature. However, at 
a given temperature, only one description is favored. By comparing the on shell actions 
for X 1 and X 2 we can obtain which geometrical description is preferred and study the 
Hawking-Page like phase transition between these two. But before we do that, we first 
study the holographic descriptions of conformal theories and demonstrate how there is no 
phase transitions. 

3.1 Conformal field theory 

For the N = 4 Yang-Mills theory on 4D Minkowski space, there are two gravity solutions 
too, one is AdS^ x S 5 (X 1 ) and the other is black hole AdS^ x S 5 (X 2 ), with the Euclidean 
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metric in Einstein frame is 



ds 



2 



r 

7? 



2 r / ^4 



n ) dr 2 + dx \ + , L ^ dr 2 + C 2 dn 2 (3.2) 
J r 2 (l- 



where is the black hole horizon and £ 4 = \kN c o! 2 where iV c is the number of D3 
branes in the dual gauge theory 11 . Note that if R,R S are the Ricci scalar in Einstein and 
string frame, we have R s = R/^fgl with = g s being the string coupling. Denoting 
string frame metric as g s MN and Einstein metric as gMN, we have Ricci scalar for the AdS 
direction as g ab R a i) ~ — Jy, a, b = 0, .., 4 and in the string frame gf'Rsab ~ —j?,L = 9s^£ 
. For supergravity to be a good approximation we need Ricci scalar to be small in string 
frame, which means we need L 2 = ^g s N c a' 3> 1. On the other hand we want to ignore 
string interactions, which means g s —> 0. Thus the limit \/g s N c a' 3> 1 can be obtained 
with N c — > oo and classical gravity with the metric (3.2) will describe the dual geometry. 



Observe that when = 0, we recover the AdS$ x S 5 metric without any black holes. 
The bulk on shell supergravity action for both X 1 and X 2 is the same 



So the action difference should come from the Gibbons-Hawking (GH) term 51, 78] 

IGH = 4- f d\^ffetl(D^) (3.4) 

where Mg is the boundary at r = 1Z, g is the metric induced on Mg, is a unit normal 
vector to Mg and is the covariant derivative. 

It is straight forward to compute the GH term for AdS§ x 5 5 with and without black 
holes. We find the GH term is 

T 2 _ 



I 2 GK = ^f- Jim (4tt*-2rjJ) (3.5) 
for black hole AdS§ x S 5 , and 

/An = km ^ 4 (3.6) 



'GH .. z 

no w 



for AdS^ x S* 5 and we have used that volume of S 5 is J dQ^ = 7r 3 . 

However, the extremal geometry X 1 only describes zero temperature. The temperature 
described by the geometry is obtained from the periodicity of Euclidean time and for X 1 
which is Poincare AdS, this period cannot be arbitrary. This is because Poincare AdS 
metric (obtained by setting = in ( |3.2[ )) is singular at r = 0. Near r = 0, with 
V\ = r 2 / C 2 the (r,r) part of the metric takes the following form 

ds 2 = — j-2 [dV 2 + V 2 V[ 2 dr 2 ) (3.7) 



"Using (2.26) with the surface integral over S 5 and F 5 = (1 + ^dh 1 A dt A dx A dy A dz, h 
readily gets N c = 
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where V( = The singularity at V\ = is really the singularity of the origin of polar 
coordinates and the metric is smooth and complete, provided [^'[(ir = d9, where is the 
polar angle. This readily gives that period f3± = J dr = 1/\V{\ f d9 — > oo as V\ — > 0. Thus 
the corresponding temperature is 



ft" 1 = (3.8) 



On the other hand, the Euclidean black hole metric is singular at r = fh- The singularity 
can be realized as that of polar coordinates with the change of variable V 2 = (l — 
and looking at the metric near the sigular point V2 = 0. This way the metric is smooth 
and complete if and only if 

h = ^~ (3.9) 

which depends on the horizon f/j. Also the temperature of the field theory at hypersurface 
p = 1Z of geometry X 2 is now redshifted due to the presence of the black hole 12 



r 1 

T 2 (K) = . p2 (3.10) 

In summary, the Poincare AdS metric without a black hole only describes the vacuum 
phase and corresponds to zero temperature. On the other hand, for any non-zero temper- 
ature, M = 4 SUSY gauge theory is described by the black hole geometry. Since there is 
no description in terms of distinct geometries for a given non-zero temperature, there is no 
Hawking-Page transition. 

Now we would like to find the thermodynamic parameters for the SU(N C ) conformal 
field theory. But since the black hole action Jq H is divergent, we need to regularize it first. 
The divergent part ~ fyTZ 4 is dependent on temperature and thus the counter term will 



also be dependent on it. Also observe that the regularization scheme adopted in [S3, 64] 
identifies AS as the regularized gravity action. Hence regularization is not just subtraction 
of the infinite part ~ fyTZ 4 , but also a finite part arising from the vacuum action. There is 



one crucial difference between the case studied in [63, 34 1 and the Poincare AdS geometry. 
Essentially, Hawking-Page considered global AdS geometry, where without a black hole, 
the geometry is regular and the periodicity of Euclidean time can be arbitrary. There, at 
a given temperature there are two descriptions, black hole and no black hole geometries 
and AS can describe the excess of energy due to presence of the black hole. For Poincare 
AdS geometry, the action S 1 only describes zero temperature, so AS is only meaningfull at 
zero temperature, where it is trivial. At any non-zero temperature, AS does not describe 
the excess energy and thus cannot be identified with the renormalized action for the black 
hole. 



12 The gtt and g cc , c = x,y,z component of the Minkowski metric at hypersurface r = 1Z are now distinct. 
The difference is dependent on 1Z and thus the local temperature at r — 1Z is dependent on 1Z 
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However, the dual conformal thermal gauge theory has free energy that scales as ~ T 4 
and entropy that scales as T 3 . To extract this information from our Poincare AdS black 
hole geometry, special care is needed while regulatizing the action. In particular, one needs 
to holographically renormalize the action following the procedure of [78] and write the 



regularized action as a functional of the induced metric at the boundary to find the stress 
tensor of the dual gauge theory. 

There is an alternative approach to obtain the thermodynamic state functions. Black 
hole geometries have entropy proportional to one quarter the surface area of the horizon 
1 79 , 8C ] . We can directly evaluate the black hole entropy for the ten dimensional AdS 5 x 



black hole geometry using Wald's formula [ 81 ]-[ 34 1 where the gravity Lagrangian is £buik 



2^2- \R — -^^r I • The result is 



27r 4 f3£V 3 1 



-AT c 2 7rV 3 T 3 (3.11) 

where T = T2(1Z = oo) and we have used definition of k\ q = 2 6 vr 7 a /4 . Knowing the entropy, 
we can of course obtain the free energy and find the appropriate counter terms that would 
lead to this free energy. The result is 

F = = -^ 3 vr 2 iV c 2 T 4 = T(/ S 2 UGRA + / 2 h + / 2 ounter ) 
T 2 _ 47^-3/2f 4 3 

^counter ~ 2 P 2 V ^ K 6 - 12 ) 

The result for entropy and free energy is consistent with the finite temperature con- 
formal field theory since they have the correct scaling with temperature. Also note that 
the entropy s = 3/4so where sq is the entropy for weakly coupled N = 4 SYM theory. We 
will not discuss this discrepancy here, please refer to [p8[ , |39[ ] and pOfl . 

Observe that the discussion about the phase transition about the conformal field theory 
is different from the case discussed in |34|]. There the field theory lives on S 1 x 5 3 , with 
circumference /3 and j3' respectively. Thus the field theory is no longer conformally invariant 
and depends on the scale P/f3'. In the large limit there can be a phase transition as a 
function of /?//?'. 

3.2 Non-conformal field theory 

Now we turn to the phase transition of non-conformal field theory. We have discussed dual 
gravity solutions of the non-conformal field theories in section 2. The zero-temperature 



solution which is discussed in section 2.1 is UV completed. The geometry is a warped 
deformed conifold in the IR and asymptotically AdS in the UV which ensures that the 
theory is holographically renormalizable. When energy flows from UV to IR as r gets 
smaller and smaller, the effective number of D3 brane is cascaded away, leaving only D5 
branes in the IR and the theory is confined. In the previous section we discussed the 
finite-temperature solution but we didn't bother to UV complete it. However notice that 
the same AdS cap can be added to the black hole geometry with a suitable intermediate 
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buffering region connecting the UV and IR region. Furthermore, far away from the black 
hole horizon, extremal and non-extremal geometry become identical, which is evident from 
the form of the metric presented in the last section. At large r S> r^, f^M < 1,1 >0 and 
thus the non-extremal metric for large r indeed becomes identical to the extremal limit. 
Hence to analyze the difference between the extremal and non-extremal geometries, we will 
ignore the effect of UV completions in this work and analyze how our results modify with 
the addition of various UV caps in our upcoming paper|K]] . 

We will now evaluate the on-shell value of the action ( |2.1[) for the extremal geometry 
X 1 without adding a UV cap and any localized sources. When the axio-dilaton field r is 
constant, the fluxes and the metric were exactly evaluated in [44|. Thus in the absence of 
D7 branes such that Nf = 0, t =constant, we get for the bulk on shell action 
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dr / dp sin#i sin#2 



(3.13) 



5g s M 2 l(p) \ 
12/i (p)Imry 

where we have used the equations of motion for G3: *ioG3 = 1C4AG3, G\ is the determinant 
of the metric ( 2.12j ) after wick rotation and G3 = — -^-H^, H^,F^ are the three form 
fluxes on deformed cone J3|. Also ho is the Klebanov-Strassler warp factor and 

pcothp — 1 



lip) 



4sinh 2 
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(sinh2 / o - 2p) 



(3.14) 



as determined in [44]. Here f3\ = T^ 1 can have any value corresponding to any temperature. 
This is because unlike the Poincare AdS geometry, the deformed cone geometry ( |2.12| ) is 
regular and after Wick rotation, the periodicity of Euclidean time j3\ can take any value. 
For a chosen value of f3\, T\ gives the local temperature of the field theory living on the 
hypersurface p = p c for any p c . 13 

Now computing the on shell bulk action for the non-extremal geometry X 2 gives 
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(3.15) 



where now the three from flux G3 is distinct from the extremal values, G2 is the determinant 
of the Euclidean metric and ph is the horizon in the radial coordinate p. By switching to 
r coordinate using the transformation (|2.23j ), the corresponding temperature of the field 
theory at hypersurface r = 1Z of geometry X 2 is 



T 2 (K) 



P~2 



02 



A-K^h{r h ) 
W(r h )\ 



(3.16) 



13 The four dimensional Minkowski metric at hypersurface p = p c has identical warp factor for both the 
time and space directions for every p c . Thus the local temperature at p = p c is independent of p c . 



- 24 - 



where is the horizon and prime denotes derivative with respect to r. The two geome- 
tries describe the same field theory at the same temperature on the hypersurface r = TZ 
if 



T = Ti(ft) = T 2 {K) =>/3i = foy/gffl 



(3.17) 



Also we have to take the GH term into consideration, this will be discussed in the next 
subsection. 

Now, to see which dual geometry is favored for a gauge theory at fixed temperature T 
and living on hypersurface 1Z = oo, we consider the action difference [E3ip3] 



AS = S\ 
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(3.18) 



where we used the relation (3.17) and <Sq H ,Sq H are the Gibbons-Hawking surface terms 
for the extremal and non-extremal actions. If AS" > 0, then extremal geometry is preferred 
and the dual gauge theory is in the confined phase. On the other hand AS < indicates 
that non-extremal geometry will be favored and the gauge theory is in the deconfined 
phase. Note that AS* only depends on p^ and thus on T. As T is altered from small to 
large values, it is possible for AS(T) to change sign. If there exists a particular T = T c 
such that AS(T C ) = 0, then T c will be the critical temperature which is realized as the 
confinement /deconfinement phase transition temperature of the dual gauge theory. Since 
we are working in the limit Nf = 0, the corresponding critical temperature is for pure glue 
large N gauge theory with no fundamental matter. 

Now to evaluate ( |3.18j ) explicitly, we need to know 5q H ,5q H ,X2 exactly, since T\ is 
already known. We will now compute *^gh whil^ analyzing X2 for large and small 
black holes separately and estimate the critical temperature. 

3.2.1 Large horizons, critical temperature and deconfined phase 

When p > is large, the internal metric simplifies from the non-extremal limit of ( 2.12j ) 
and takes the form ( p. 32 ) where the radial coordinate r given by ( p.23 ) becomes convenient. 
Thus for large radial distances we will use r coordinate and we really have a Klebanov- 
Tseytlin type geometry. Then the three form fluxes in the presence of a black hole are 



14 The temperature (3.16) is derived by taking the non-extremal metric to be (2.32) with a(r) = l. That 
is we are assuming the horizon is large enough such that the non-extremal limit of (2.13) takes the form 
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given by (gjg),(HD while the flux strength is given by (|2.55[ ), Furthermore, only keeping 
terms up to linear order in 0(g s M 2 /N) and using = g s , I2 can be exactly evaluated 
with the following result 
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On the other hand for large p, we can use the r coordinate to get 

5r 3 g s M 2 sin#i sin#2 
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(3.20) 



To compute the GH term for KS geometry with black hole, we find the embedded Euclidean 
metric at the boundary, 
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dsl 
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+ o/Jn(i + *^^]ogL + Ai(r)) (b(r)g$dxW) 
where A\ (r), G(r) has the form given by ( 2.64[ ) . The GH term is 
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where g is the metric ( |3.21| ). Since the internal metric g pq = k{r)g P q is unknown, we 
can not determine the GH term completely. However, from eq. ( |2.64| ) we know that the 
corrections to g pq is of order 0(g s M /N)-^,l > 0, so their contributions to the GH term 
will be of same order. Using this scaling, we find that the GH term for KS black hole 
metric is 
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where Vg = V3 x V5. Also the term ~ d r^gsM 2 /N comes from the g' pq ,Ai(r) and G(r) 
corrections that modifies the unit normal vector n M and its covariant derivative. In deriving 



(3.23), we only kept terms up to linear order in 

log f^V When there is no black hole, i.e rh 
geometry. That is 
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which is valid when iV 3> M, N S> 
= 0, we get the GH term for KS 
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(3.24) 
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For large horizons, p > ph is always large and thus we can use ( |3.20 ) near the horizon 
to obtain the action difference 

AS = °>*W*j Hm A ( K\ _ 59 _ 27( , _ al) \ (3 25) 

where a 1 arises from the following expansion 

?r , x #fc ( 729a 1 g s M 2 \ f ? „ „ f g s M 2 \ 

e 2B = g(r) = 1 - ^ (i + ^ j + _| ^«L_ j + ... ( 3 .26) 

Due to the presence of D5 branes, we expect the black hole horizon to be larger than 
the Schwarzschild value |57|| , which is possible if a 1 > 0. In the rest of the analysis, we will 

assume a 1 > 0. Note that d is determined by the Einstein equations and the flux equations 

f 4 

once boundary conditions are imposed. Since d arises from the coefficients of the term -| 
in the expansion ( 2.64| ) of g' pq , Ai(r), G(r), it is sensitive to the near horizon geometry and 



in particular can be obtained from the horizon values of these functions. Furthermore, d 
is related to Ai(r/j), which in turn determines the number of effective degrees of freedom 
at a temperature T ~ r^/L 2 . Thus from the gauge theory side, d is related to the effective 
colors at the thermal scale. As we are not able to exactly solve all the Einstein equations 
near the black hole horizon, we can speculate two scenarios: 

• d < -(59/108 - a 1 ): In this case AS > 0, which means extremal geometry is favored 
and black hole is not formed. In that case all black holes have higher free energy and for 
all corresponding temperatures, the gauge theory will be in the confined phase. However, 
this does not mean there is no phase transitions at all- rather upto linear order in 9s ^ , 
AS > 0. In fact there are terms in AS which are second or higher order in 9s ^ log (jf^J 

and for TZ — > oo, it is possible that log ( — ) 3> 1. Then the perturbative analysis 



breaks down and we cannot ignore higher order terms. When the higher order terms are 
accounted for, A5 could indeed be zero for some r£ and there would be a phase transition. 

• d > -(59/108 - a 1 ): In this case, using ( |3.25| ), it is possible to obtain AS" = with the 
following value for critical horizon 

TZ 

f% = exp(f + 27(d-ai)) ^ 

There is always a phase transition provided 9s j^ log (j^J ^ 1- This is possible with the 
scaling TZ = N l l i ^/~a' — > oo, r\ ~ 0(\/a7) and N S> M 2 . Then the corresponding critical 
temperature is 
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where A = g s N is the 'tHooft coupling. For T > T c , the black hole is formed and describes 
the deconfined phase which has non-zero entropy. On the other hand for T < T c , we 
have a confined phase with zero entropy described by the extremal geometry. Thus at 
T = T c , there is a first order phase transition in the gauge theory and we have obtained 
a gravitational description of it in terms of Hawking-Page phase transition between two 
geometries. 

Now to obtain thermodynamic state functions of the gauge theory, we first need to 
renormalize the action. From eq. ( 3.1S| ) and eq. ( |3,24 ), we see that the on shell action is 



divergent and we need the following counter terms: 

r2 ^4 



counter ^ y 2J ' U8N 1 54 ft4 J 



(3.29) 



where k can be determined by matching the Bekenstein-Hawking entropy of the black hole 
to the entropy obtained through the partition function using the relations (||1| and (|2~Tll ). 



With this counter term, for T > T c , we can write down the free energy and entropy for the 
deconfined phase of non-conformal field theory using ( [2.10 , |2,11 ) 
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~ -iV 2 T 4 (l + ^log(LT) 

s ~ N 2 T 3 (l + (1 + log(LT))) (3.30) 

where b is some constant independent of M, N but determined by the relation between 
fh and T. We observe that —F/rt > for N 3> M and increases with r^. This indicates 
entropy increases with tv Also temperature T is given by ( 3.1 6| ) and since the exact metric 



at the horizon is not known, the exact expression for T as a function of fh, M and N is also 
unknown. This is the reason the exact scaling of free energy with temperature is unknown 



and in (fi.30|), we approximated T ~ — *%n, which is a reasonable since it is the leading 
term and g s M 2 /N -C 1. Using the same renormalization procedure for the extremal action 
gives -^counter = ^counter (A — ^ /^2 5 fh = 0) an d we get F = 0, s = for T < T c - the confined 
phase. 

The above analysis, done purely using supergravity can be interpreted as a possi- 
ble proof of deconfinement at large temperatures for a gauge theory which has a gravity 
dual. On the other hand, when temperature is zero we have the unique extremal geometry 
corresponding to the confined phase. Thus starting with a confining phase at zero temper- 
ature, our analysis shows that the gauge theory transitions to a deconfined phase provided 
d > —(59/108 — a 1 ) i.e the deformations of T ' near the horizon behave in a particular 
way. Thus, we have a confinement/deconfinement phase transition and the entire analysis 
is done using the low energy limit of type IIB superstring action ([2.1]), 

3.2.2 Small horizons and the confined phase 

From ( 3.25D , observe that as long as < rf, A 5 > and the extremal geometry without 



a black hole will be preferred over the black holes geometry. Thus when the horizon radius 
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is small compared to the critical radius the dual gauge theory will be in the confined 
phase. 

Note that our entire analysis has been done in the large p region where the relevant 
radial coordinate is r. Going back to the p coordinates, we see that the critical horizon in 
terms of p coordinate becomes 

Pc ~ 31og (3.31) 

The horizon p^ ~ 31og < Pc can still be a large number referring to a large black hole. 

But even this 'large' black hole is not favored against the extremal solution since pt < p c - 
Thus critical horizon is the relevant scale that determines largeness and smallness. 

For ph <C Pc, we do not have explicit form of the metric and the flux and thus cannot 
evaluate the on shell non-extremal action S 2 . Equation ( |3.25| ) does not apply in such cases 
and we need explicit solutions of the metric and fluxes to evaluate the action difference. 
Note that for these small black holes, we have to consider the non-extremal limit of ( 2.12| ), 



which is not known. The difficulty arises for the following reason: the proposed solution 
Q2.32 ) with the scalings ( |2.64 ) is only valid for large black holes. In finding these scalings, 



we heavily relied on the smallness of the parameter g s M 2 /N. For small black holes, such 
scalings do not apply since at small p, N cascades down to smaller and smaller values. 
The relevant quantity now is g s M 2 /N eS (p) which is no longer small and the perturbative 
analysis can no longer be applied. 

However, using ( [3.25 ) which is only applicable for large i.e large ph, we observe that 



AS = only once. On the other hand small p^ means small temperature and in particular 
at zero temperature, the SU(M) pure glue dual gauge theory confines. Since there is 
nothing special that occurs in the gauge theory at low temperature, it is reasonable that 
AS 7^ at low temperature. This implies AS becomes zero only at high temperatures and 
we have unique critical temperature given by the critical horizon p c in (3.31). For every 



Ph < Pc, no matter how small it is, we will have confinement. Of course this is consistent 
with confinement at zero temperature and the scenario that deconfinement happens only at 
a single large temperature- indicating we have only two phases, confined and the deconfined 
phase. To gain a better understanding of the gauge theory and its precise degrees of freedom 
in the different phases, we will now sketch a possible brane configuration that gives rise to 
the dual geometries X 1 and X 2 . 



4. Brane configurations, the gauge theory and connection to QCD 

In the previous sections we studied two gravity solutions and the Hawking-Page like phase 
transition between them, in this section we will study the brane configurations that lead 
to the dual gauge theories. 

Suppose we place N D3 branes at the tip of a regular cone [See Fig 1(a)]. At zero 
temperature, the gauge group is SU(N) x SU(N) with bi- fundamental fields Ai,Bj,i,j = 
1, 2. This gauge theory has a conformal fixed plane and the number of D3 branes remains 
the same at all energy scales. Now if we put another stack of D5 branes that wraps 
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the vanishing two cycle at the tip of the cone [See Fig 1(b)] the gauge theory becomes 
SU(N + M) x SU (N) with the bi-fundamental fields, and it is no longer conformal. The 
SU(M + N) sector has 2N effective flavors while the SU(N) sector has 2(N + M) effective 
flavors thus it is dual to the SU (N — M) x SU (N) gauge theory under an Seiberg duality. 
Under a series of such dualities which is called cascading, at the far IR region the gauge 
theory can be described by SU(M) x SU(K ) group, where N = IM + K , I, < K < M 
are positive integers. Now the number of 'actual' D3 branes N is no longer the relevant 
quantity, rather N ± pM where p is an integer describes the D3 brane charge. We take 
K = in all our analysis, so at the bottom of the cascade, we are left with M = 1 SUSY 
SU(M) strongly coupled gauge theory which looks very much like strongly coupled SUSY 
QCD. 




(a) N D3 branes placed at the tip(b) M D5 branes wrapping the van-(c) M Anti-5 branes separated from 
of the conifold ishing two cycle of the conifold the D5 and D3 branes but also lo- 

cated at the tip r = 0. 



Figure 1: Brane construction of conformal field theory, non-conformal field theory and UV com- 
pleted non-conformal field theory 



Due to the strong coupling at the IR, the gauge theory develops non-perturbative 
superpotential j73| and breaks the Z<im R symmetry down to Z2 group. Since the complex 
fields A{,Bj,i,j = 1,2 also describe the complex coordinates of the cone, the breaking of 
the R symmetry modifies the geometry from a regular to a deformed cone [|44|j . Thus to 
capture the IR modification of the gauge theory, we must consider the warped deformed 
cone. 

All the discussions so far have been concerned with zero temperature dynamics of 
the gauge theory. At finite temperature, we expect the number of D5 branes M remains 
constant at different temperatures, thus on the gravity side the units of F3 fluxes remains 
the same too. The three form flux F3 is given by ( 2.5C| ). This is the crucial step in our 
analysis since in the presence of the black hole, the base T ' is modified and the spheres 
S 2 ,S 3 are now squashed and deformed. But we assume that at r 1— > 00 the base is still T ' 
and the S 3 at the infinity still encloses the directions which the D5 branes extend along on 
the gauge theory side. Thus the integral of the RR flux F3 over this S 3 should be equal to 
the the number of D5 branes M and F3 is the same as in the extremal case. 

Although F3 remains the same with or without the black hole 15 , #3 given by ( 2.53j ) up 
to linear order in 0(g s M 2 /N) is distinct from the extremal flux, H® = 3ga ^/" dr AW2. So 

15 -F3 remains the same but F3 = F3 — C0H3 changes in the presence of the black hole, which is consistent 
with our earlier analysis J57{ |. 
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the three form flux G3 in the presence of the black hole is modified and it is no longer ISD. 
Since #3 is modified, B2 is also modified and this means that the change in B2 can alter the 
RG flow. But as discussed in 57 j, this modification does not alter how the gauge theory 
couplings flow with energy scale, it only redefines the energy scale in terms of temperature. 

Now observe that the gauge group becomes SU(M) only at the far IR while at high 
energies, it can be described by SU(k(A)M) x SU((k(A) — 1)M) group, with k(A) in- 
creasing with energy. Thus the UV of the gauge theory has divergent effective degrees of 
freedom and looks nothing like QCD. Although the confined phase of the gauge theory 
may resembles ftf = 1 SUSY QCD, the deconfined phase of the gauge theory is quite dif- 
ferent. Hence the value for critical temperature ( |3.28 ) cannot be directly related to the 
confinement /deconfinement temperature of QCD. 

As already discussed in section [O], to make connections to large N QCD, we must add 
localized sources that can modify the UV dynamics of the gauge theory and consequently 
alter the large r region of the dual geometry. One possibility is to add M anti-five branes 
separated from each other and from the D3 — D5 branes at the tip of the cone. To obtain 
this separation, we must blow up one of the S 2, s at the tip and give it a finite size - which 
essentially means putting a resolution parameter. The setup is sketched in Fig. 1(c) [we 
did not blow up the tip, which was intentional and will be explained in what follows]. The 
separation gives masses Ao to the D5 — D5 strings and at scales less than the mass, the 
gauge group is SU(N + M) x SU(N) x U(1) M where the additional U(l) groups arise 
due to the massless strings ending on the same D5 brane. At scales much larger than Ao, 
D5 — D5 strings are excited and we have SU (N + M) x SU(N + M) gauge theory. For 
A < Ao, i.e. at low energy, gauge theory is best understood as arising from the set up of 
Fig. 1(b) (since the modes from D5 branes are not excited). At high energy A 3> Ao, the 
gauge theory is best described as arising from Fig. 1(c). Essentially at high energies, the 
resolution ro ~ Ao is not significant and we have effectively N + M D3 charge at the tip 
of the regular cone, which arises by placing M number of D5 — D5 pairs at the tip. 

Of course there will be tachyonic modes arising from D5 — D5 strings and one needs 
to consider electric and magnetic fluxes on the respective branes to stabilize the system 
|85]-[91|. Furthermore, to incorporate fundamental matter and chiral symmetry breaking, 
we need to add D7 — D7 branes to the brane configuration of Fig. 1 and the details were 
discussed in 58 ]. The RG flow arising from a UV complete model with D3,D5 and anti 
D5 branes in sketched in the Fig. 2. 

On the supergravity side, the anti five branes will source RR and NS-NS three form 



flux such that the total flux F^,H^ — > as r — > 00. As already discussed in section 2.1 



vanishing of the fluxes indicate that the gauge theory reaches conformal fixed point where 
the two gauge couplings 51,52 becomes identical. At the far UV, since both the groups 
have same rank, we effectively have a single Yang Mills coupling ^ym — > with 'tHooft 
coupling <?YMAeff held fixed at large value. Thus this UV completion arising from anti 
brane sources in principle gives rise to a QCD like theory which confines in the IR and 
becomes conformal in the UV. In fact the Yang Mills coupling of the gauge theory becomes 
free in the far UV- just like QCD, indicating that the gauge theory is indeed very similar 
to large N QCD. 
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SU(N+M) XSU(N+M) 
g r 1/A" 

Asymptotic AdS geometry 



; Flow from UV to IR 

f 

SU(N+M) X SU(N) 
g~log A 

Klebanov-Strassler type geometry 



; Flow from UV to IR 

! 

SU(M) 
gr-log A 

Klebanov-Strassler type geometry 



Figure 2: Gauge theory arising from the brane configuration in Fig. 1(c) 
5. Conclusions 

In this paper we studied the non-extremal geometries on the dual gravity side of the 
non-conformal finite temperature field theory. To find the exact non-extremal geometry 
analytically is extremely difficult because the solution is no longer supersymmetric and we 
have no control over the internal manifold. Even numerically it will be a formidable task 
to solve all the Einstein equations and field equations together. However up to linear order 
in 0(g s M 2 /N), the supergravity equations drastically simplify and the deformations of the 
internal metric along with the corrections to black hole factor and the warp factor can be 
evaluated as an infinite Taylor series in Using this expansion, it is straight forward 
to evaluate the on shell ten dimensional gravity action and we find a phase transition 
between the extremal and non-extremal geometry. Since the geometry is dual to a non- 
conformal gauge theory, this transition is interpreted as the confinement /deconfinement 
phase transition. 

The non-extremal gravity solution we discussed is not UV completed. However, under 
the scaling 1Z = N 1 ^ 4 \/ r a' — > oo and N 3> M, we get 9s ^f log^ to be small and there is 
no logarithmic divergence. In fact there are no log(7£) divergence in the Gibbons-Hawking 
boundary term. The only term that has logarithmic term arises from the bulk action- but 
since fh > scales with 7Z, the term log^ is in fact convergent. The critical temperature 
of the confinement/deconfinement we find in section || is obtained using this particular 
scaling and thus our analysis is somewhat restrictive. It will be interesting to UV complete 
the geometry and then there will be no need for this particular scaling since there will no 



UV 



Scale A 



IR 
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logarithmic divergences to begin with. Using the UV complete geometry to study the phase 
transition and other thermal properties of the non-conformal field theory is the content of 
our upcoming work 6C]. 

We also did not include effect of running dilaton in the black hole geometry and 
did not consider D7 branes and other localized sources in non-extremal dual geometry. 
This means the thermodynamics we obtain is strictly restrictive to gauge theory with no 
flavor, no Baryochemical potential and it is not surprising that we obtain a first order 
phase transition. However the perturbative procedure outlined here to solve the Einstein 
equations along with the flux equations can easily be generalized in the presence of a 
running dilaton field and other localized sources. Similarly the on shell action in the 
presence of running axio-dilaton fields and localized sources can also be evaluated up to 
linear order in our perturbative parameter. Note that for a UV complete scenario, we 
must replace M — > M (r ), lim^oo M(r) — > and thus the on shell value of the action 
will be significantly different. In fact, the presence of localized sources will dictate how 
fast M(r) vanishes- which also indicates how rapidly the theory becomes conformal. This 
means the width of the conformal anomaly will be highly sensitive to the details of the 
localized sources. In [58] a detailed gauge theory analysis was done and the non-extremal 



dual geometry will be presented in [6C] 
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A. Appendix 

In this appendix we study a = e iA+B case in section 2, assuming the internal 5D manifold 
M is T 1 ' 1 and G 3 = 0. 

By simplifying the Einstein equations and flux equations we finally get the following 
four equations, 

k 5 e- 4A+B (4A' + B') = 1 

2B " + AB' 2 + 5^ - ^B' - + 10 jB' = 

2e 2B akk" + 8e 2B ak' 2 + Ae 2B akk'B' + Ae 2B ak 2 B' 2 

-e 2B a'k 2 B' - e 2B a'kk' - 8a 2 + 2e 2B ak 2 B" = 



a . _, b" 



a 10 ¥- 4B '- 2 1F= (A1 ' 



After further simplification we find, 



1 , f HsL 

k= (l- e 2B)l/4 ' a = e J fc ' 

B' + 2 - 2e- 2B = 0, or B' 2 {1 + Ae~ 2B ) = 2B"(1 - e~ 2B ) (A.2) 
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The first differential equation of B gives e 2B = 1 — 1/r 4 which is the usual AdS black hole 
solution. The second differential equation of B can be simplified as 

B' = e B (-l + e- 2B ) 5 /2 (A.3) 

It does not have an analytic expression, but we can see that e 2B < 1. When B <C eq. 
( |A.3| ) can be approximately solved e 2B ~ 4r. Once B is known we can easily get A, a and 
k. 

Now if G3 is non-vanishing it should be ISD as in the GKP paper, and thus it does 
not effect the internal manifold at all, the only thing that changes is the Bianchi identity 
of the five form flux, and thus the warp factor A. So we have another solution with the 
same o, k and B but different A. 

We are not very sure about the use of this kind of solutions or whether they are gravity 
duals to some field theory. It might be interesting to get numerical solutions and explore 
the properties of them such as stability, KK reductions, etc. 
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